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Abstract. Bethe-Salpeter (BS) equation in Minkowski space for scalar particles is solved for a kernel given 
by a sum of ladder and cross-ladder exchanges. The solution of corresponding Light-Front (LF) equation, 
where we add the time-ordered stretched boxes, is also obtained. Cross-ladder contributions are found to 
be very large and attractive, whereas the influence of stretched boxes is negligible. Both approaches - BS 
and LF - give very close results. 
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1 Introduction 

In a preceding paper we have proposed a new method 
for solving the BS equation [2] in the Minkowski space. 
Our approach does not make use of the transform to the 
Euclidean space (Wick rotation) and is applicable to an 
arbitrary interaction kernel. This method was first applied 
in pP to obtain the ladder solutions of a scalar model 
with interaction Lagrangian Cwc — 9<t> X- F° r niassless 
exchange we reproduced analytically the Wick-Cutkosky 
equation 0. For massive ladder exchange our numerical 
results coincide with ones found in previous works based 
on the Wick rotation. 

In the present paper we solve the BS and LF equa- 
tions with a kernel given by a sum of ladder and cross- 
ladder graphs. This constitutes the first calculation of 
cross-ladder effects in both equations. We consider scalar 
constituents in the state with zero total angular momen- 
tum. BS equation with ladder kernel was solved in Min- 
kowski space in 0]. 

Non ladder effects, within the same model, using Feyn- 
man-Schwinger representation, were considered in ref. [S]- 
In this work the full set of all irreducible cross-ladder 
graphs in a bound state calculation was included. In jS] the 
effect of the cross-ladder graphs in the BS framework was 
estimated with the kernel represented through a disper- 
sion relation. Non ladder self-energy effects in LF equation 
have been incorporated in [7||H]- 

The plan of the paper is the following. In section [21 we 
sketch the method used in for solving the BS equa- 
tion and find the kernel corresponding to the cross-ladder 
Feynman amplitude. In section [3] we remind the LF equa- 
tion, obtain the LF cross-ladder kernel and incorporate 
also two stretched box graphs. Numerical results are given 
in section^ Section [S] contains some concluding remarks. 
Technical details of the derivation are given in appendices 
ElElandO 



2 Cross-ladder kernel in BS equation 

The method for solving the BS equation proposed in our 
previous work L is based on projecting the original equa- 
tion: 



<£(£;, p) 



[(§ + fc) 2 - m 2 + ie] [(§ -fc) 2 
■iK(k,k',p)<P(k',p) 



d 4 k' 



(i) 



on the LF plane, i.e in applying to both sides of the 
integral transform, which for the left-hand side reads 



$(k + [3uj,p)dP 



and similarly for the full right-hand side. Here w is an 
arbitrary four- vector with u 2 = 0. This transformation 
eliminates the singularities of the BS equation. The BS 
amplitude in the transformed equation is then written in 
terms of the Nakanishi integral representation jS], which 
for zero angular momentum reads: 



$(k,p) 



dz 



d'y 



(2) 



[7' + m 2 



\M 2 



k 2 — p ■ k z' 



The equation satisfied by the weight function 3(7, z) was 
derived in ,1, and has the form: 



g{Y,z)dy 



7' + 7 + z 2 m 2 + (1 — z 2 )k 2 



dj' / dz'V( 7 ,z; 7 ',z')g^',z'), 



(3) 
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where V is a kernel given in terms of the BS interaction 
kernel iK by 



v(rr,zn',z') 



oj ■ p 



-il(k + /3uj,p)dl3 



[(! + *• 

1 



m 2 + ie\ 



[(§ - k - (3ujf - m 2 + ie] ' 
d A k' iK(k,k',p) 



(4) 
(5) 



( 2n ) [k' 2 + p ■ k'z' - i - k 2 + ie] 
The bound state mass M enters through the parameter 

k 2 = m 2 - -M 2 . 
4 

Equation J3J is equivalent to the initial BS equation 
and provides, for a given kernel K(k,k',p), the same 
bound state mass M . Once g(-y, z) is known, the BS am- 
plitude can be restored by eq. Corresponding LF wave 
function ip(kj_,x) can be easily obtained by 



x(l - 1)5(7' , 1 ~ 2x)d7' 



V^->0 7 ' + fc 2 + m 2 _ x (l _ X ) M 2 

In P for the ladder exchange 

„2 

X (i) (fc,/c',p) 



(6) 



A;') 2 — (i 2 + ie 



we calculated the integrals (@J, (J5J for the kernel V and 
solved equation ©. Derivation of V for the cross-ladder 
kernel is quite similar but more lengthy, since the kernel 
itself is more complicated. The cross-ladder BS kernel is 



Pi 



In the bound state problem the value of p 2 is fixed: p 2 = 
M 2 . 

The amplitude corresponding to the diagram in fig. ^ 
reads: 

K [CL \k,k\p) = 

^4 



(27r) 4 J [p" 2 - m 2 + ie] \{p' 2 -pi+ p") 2 - m 2 + ie 



4„// 



d 4 p 



[(Pi ~ P") 2 - V 2 + ie] [01 - P") 2 - M 2 + « 



(7) 



We have first to calculate this expression, substitute the 
result in (J5J, then in Q and find in this way the cross- 
ladder contribution to the kernel V(7, 2; 7', z') in equation 
©■ 

The full kernel - including ladder and cross-ladder 
graphs - will be written in the form: 

V (7, 2; 7', *0 = (7, 2; 7', A + V^ CL) (7, ^; 7', A 

The ladder kernel was found in pp. The cross-ladder 
contribution V^ CL ' is calculated in appendix El 



3 Ladder, cross-ladder and stretched-box 
kernel in LF equation 

We would like to compare the results obtained in the BS 
approach with the equivalent ones found in Light-Front 
Dynamics (LFD). For this aim we precise in what follows 
the LF equation and derive the corresponding kernel. In 
the well-known variables fcj_ and x the LF equation reads 
(see e.g. [TU]): 



a;(l — x 
2 



(8) 



p^p 



m 2 f , , is , , , 9> d k\ dx' 

J i>(k> ± ,x>)V LF (k' ± ,x';k ± ,x,M 2 ) 2xl ^_ xl) 



p 2 -p, + p 
Fig. 1. Feynman cross graph. 

shown in fig.^ On mass-shell it depends on two variables: 
s = {pi +P2) 2 and t = (pi -pi) 2 . The kernel K in the BS 
equation |QJ is off-mass shell. It depends also onp 2 , i.e., on 
six scalar variables in general: s, i,p 2 ,Pa,pi 2 \p'^ ■ One can 
also construct six variables, using the total momentum p 
and relative momenta k, k': 



k a 



where 

P=Pl+P2 =pi+P2, k 



k ■ k , k ■ p, k ■ p, 



■(P1-P2), k' = -01 -Pa)- 





a 

k ' z K 




k, 


k', k, 


k', 









Fig. 2. Ladder LFD graphs. 



The two time-ordered ladder graphs are shown in fig. 
[2] The corresponding kernel has the form: 



V&\k' ± ,x';k ± ,x,M 2 ) = 

A-na9(x' - x) A-kolQ(x - x') 

~ 'Jx> - x){s a - M 2 ) ~ (x — x')(st — M 2 ) 



(9) 
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where 



k\ + m 2 (k'± - k\_) 2 + /i 2 k'\ + m 2 
x x' — x 1 — x' 



and 



k' 2 , + m 2 (k'j_ - kj_) 2 + a 2 k 2 + m 2 



1-x 



The six different cross-ladder diagrams are shown in 
u g- El They have order a 2 . In addition, and in the same or- 
der a 2 , there are two time-ordered irreducible graphs with 
two mesons in the intermediate state (stretched boxes). 
They are shown in fig. 01 The full LFD kernel - includ- 
ing ladder, cross-ladder and stretched-box graphs - will 
be written in the form: 



V L F(k' ± ,x'; k±,x, M 2 ) = V [ LF ] + V^> + 
The term V^p is given by © and 



r(CL) 



r (SB) 



(CL) 



LF 



(10) 



i=l 



i=7,f 



are calculated in appendices iBl and [HI correspondingly. 



4 Numerical results 

We have solved numerically the BS equation in the form 
J3Jl and the LF equation (JHJ both for the ladder (L) and 
ladder+cross-ladder (L+CL) kernels. In the case of the LF 
equation (JSJ we added also the stretched box contributions 
(L+CL+SB) with two intermediate mesons shown in fig. 
01 These stretched box contributions, as well as those with 
any number of intermediate mesons, are generated by it- 
erations of the ladder BS kernel and are thus implicitly 
included in the BS approach. 

The numerical procedure, based in the spline expan- 
sion of the solution, is similar to the one used in pp. By 
expanding the solution g(j, z) on a spline basis we obtain 
a matrix equation: 



XB(M)g = A(M,a)g. 



(11) 



Like in pQ, the matrix B(M) was regularized by adding 
to its diagonal a small value e ~ 10~ 4 1CP 12 and we 
have checked the stability of the eigenvalus relative to e. 
The difference with respect to the ladder kernel is that the 
coupling constant a does not appear linearly in the right- 
hand side of For a fixed mass M, the eigenvalue A 
is calculated for different values of a and the value cor- 
responding to A = 1 can be easily extrapolated from the 
almost linear behaviour of A (a). 

The binding energy B as a function of the coupling 
constant a is shown in figures El and El for exchange masses 
/i = 0.15 and fi — 0.5 respectivley and unit constituent 
mass (m = 1). Corresponding numerical values are given 



in tabled Obtaining this results is quite a lengthy numeri- 
cal task, entirely due to the evaluation of the 4-dimensional 
integral required in the CL kernel. This fact makes diffi- 
cult - though straightforward - to reach the same accuracy 
than for the ladder case. The results in tabled have been 
obtained with N as 10 gauss quadrature points in each 
integration variable ensuring an accuracy of 1%. 

We see that for the same kernel - ladder or (ladder + 
cross-ladder) - and exchange mass - [i = 0.15 or /i = 0.5 - 
the binding energies obtained by BS and LFD approaches 
are very close to each other. The BS equation is slightly 
more attractive than LFD. At the same time, the results 
for ladder and (ladder +cross-ladder) kernels considerably 
differ from each other. The effect of the cross-ladder is 
strongly attractive. Though the stretched box graphs are 
included, its contribution to the binding energy is smaller 
than 2% and attractive as well. This agrees with the direct 
calculation of the stretched box contribution to the kernel 
done in 

The zero binding limit of fig. El deserves some com- 
ments. It was found in J21 that for massive exchange, 
the relativistic (BS and LF) ladder results do not coincide 
with those provided by the Schrodinger equation and the 
corresponding non relativistic kernel (Yukawa potential) 
even at very small binding energies. Their differences in- 
crease with the exchanged mass /j, and do not vanish in 
the limit B — > 0. We have displayed in fig. [7| a zoom of 
fig. El for small values of B. We see from these results that 
the cross ladder and stretched box diagrams reduce the 
differences but are not enough to cancel it. 



— • LF Ladder 

- - - LF Ladder+CL+SB 

Schrodinger (Yukawa) 




0.03 



Fig. 7. Zoom of figure |S| in the zero binding energy region. 
The ladder and (ladder +cross ladder +stretched box) results 
obtained with the LF equation are compared to the non rela- 
tivistic ones (Schrodinger equation with Yukawa potential). 
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Fig. 3. Cross LFD graphs. 
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Fig. 4. Stretched boxes. 
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Fig. 6. The same as in fig. for exchange mass /i = 0.5 and, in addition, binding energy B for LF equation with the ladder 
+cross-ladder +stretched box (L+CL+SB) kernel. 
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Table 1. Coupling constant a for given values of the binding energy B and exchanged mass n — 0.5 calculated with BS 
and LF equations for the ladder (L), ladder -fcross-ladder (L+CL) and (in LFD) for the ladder -(-cross-ladder +stretched-box 
(L+CL+SB) kernels. 



B 


BS(L) 


BS (L+CL) 


LF (L) 


LFD (L+CL) 


LFD (L+CL+SB) 


0.01 


1.44 


1.21 


1.46 


1.23 


1.21 


0.05 


2.01 


1.62 


2.06 


1.65 


1.62 


0.10 


2.50 


1.93 


2.57 


2.01 


1.97 


0.20 


3.25 


2.42 


3.37 


2.53 


2.47 


0.50 


4.90 


3.47 


5.16 


3.67 


3.61 


1.00 


6.71 


4.56 


7.17 


4.97 


4.91 



5 Conclusion 

We have solved, for the first time, the BS equation for the 
kernel given by sum of ladder and cross-ladder exchanges. 
The solution was found in Minkowski space, i.e. without 
making use of the Wick rotation, by a new method devel- 
oped in pp. 

In order to compare two different relativistic approa- 
ches, we have also solved the corresponding LF equation. 

We have found that the cross-ladder contribution, rel- 
ative to the ladder one, results in a strong attractive ef- 
fect. The BS and LFD approaches give very close results 
for any kernel, with BS equation being always more at- 
tractive. These approaches differ from each other by the 
stretched-box diagrams with higher numbers of interme- 
diate mesons. Our results indicate that the higher order 
stretched box contributions are small. This agrees with 
direct calculations in LFD of stretched box kernel (fig. 
0} with two-meson states and with calculations of 
the higher Fock sector contributions J2| m the Wick- 
Cutkosky model. Calculation in LFD of binding energy 
with the stretched box contribution (L+CL+SB) and its 
comparison with (L+CL) also shows that the stretched 
box contribution is attractive but small. 

The comparison of our results with those obtained in 
jS], evaluating the binding energy B a u for the complete 
set of all irreducible diagrams, shows that the effect of the 
considered cross ladder graphs, though being very impor- 
tant, represent only a small part of the total correction. 
Thus for \i — 0.15 and a = 0.9 the corresponding bind- 
ing energies obtained with BS equation are Bl 0.035, 
Bl+cl » 0.06 and B aU « 0.225. 

We are aware about the only paper 6 where the sep- 
arated effect of the cross-ladder in the BS framework has 
been estimated. The method is based on writing an ap- 
proximate dispersion relation for the kernel. Our results 
are smaller than the ones found in this reference by a fac- 
tor 3. 
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A Calculating BS cross- ladder 

We calculate here the cross-ladder contribution V^ CL ^ to 
the kernel in eq. J3J| . We start with the cross-ladder am- 
plitude, eq. J7|). Using the Feynman parametrization: 

i /-i ri—y3 r^—yi—vz 

6 / dy 3 / dy 2 / dyi 
Jo Jo Jo 

1 



abed 



[ayi + by 2 + cy 3 + d(l - yi - y 2 ~ J/3)] 4 

and calculating then the integral (0) over p" in a standard 
way by means of the formula 



d A p 



{p" 2 + A + ie) n (n - l)(n - 2) (A + ie) n ~- 
we find: 

K^ CL \k,k',p) = 



^ (12) 



(13) 



16a 2 m 4 / dy 3 / dy 2 



"1—3/3 pl—V2—V3 ^ 

dyi 4^ 

Jo Jo Jo A 

where a — g 2 / (I6irm 2 ) and 

A' = k' 2 (yi + y 3 )(l - y 1 - y 3 ) 
+ k 2 (y 2 +y 3 )(l ~y 2 - y 3 ) 

+ k' -p (yi(l - 2/1 - 2/2) - 2/3(2/1 + 2/2)) 
-k-p (2/2(1 - 2/1 - 2/2) - 2/3(2/1 + 2/2)) 
+ 2k-k' ((1/1 + 2/3) (2/2 + 2/3) ~ 2/3) 

+ Jm 2 (2/1+2/2X1-2/1-2/2) 

- P 2 (1 - 2/1 - 2/2) - m 2 (2/1 + 2/2) + it- 



(14) 
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Substituting the kernel in eq. (JjjJ, we obtain: 
d A k' iK (CL ){k,k',p) 



1 rl— V3 fl— yi— V3 

dV3 / dy 2 
o Jo Jo 



(15) 



dyi 



d 4 k' 16ia 2 m 4 



(2tt) 4 A' 2 m ' 
where A' is given by (|14fl and 77 is denominator in JSJ: 

1 



H = k + p ■ k! z' 

Using the formula 
1 



~-M l 



7 + le. 



A' z H 3 



12y 4 (l - Vifdyi 



y 4 A' + (1 - Vi)H 



and shifting the integration variable k! to eliminate the 
terms linear in k', we calculate, again by means of eq. 
I(12|). the integral (|T5|) over fc'. The result is represented in 
the form: 



I(k,p) = - 



a 2 m 4 f riy^l ~ y i ) 2 dy 1 dy 2 dy 3 dy i 



A(k,p) + ie 



where 



V = 1 - J/4[l - (1 - Vi ~ 2/3)(?/i + 2/3)] 



(16) 



(17) 



and A(k,p) depends on k,p, 7', z', 2/1-4- We separated the 
factor 77 in numerator of (|TBj) so that A(fc,p) be a polyno- 
mial in all the variables. We do not precise it. 

Now we shift the argument k of I(k,p): k — > k + (3lo, 
substitute I(k + (3w,p) in eq. Q for V and replace p" = 
(to ■ p)f3. In addition to the variables 7', z', which enter 
through eq. JSJ, the kernel V depends on three scalars k 2 , 
p-k and ^:V = V{k 2 ,p-k,^;i,z'). They vary in the 
intervals 

- 00 < k 2 < 0, 



1 uj-k 1 
-00 < p ■ k < 00, < < — 

2 ~ uj-p ~ 2 



and satisfy the relation 



p- fc = 2 (fc 2 

U! ■ P 



(18) 



(19) 



Therefore, only two of them are independent. We intro- 
duce two new variables 7, z related to the three old ones 
as: 



fc 2 

p ■ k 

lu ■ k 
to ■ p 



(7 + z 2 m 2 ) 
1-z 2 ' 
z[y + z 2 m 2 + (1 — z 2 )k 2 ] 



1 



1 

'2 Z - 



With these definitions, the relation (|19f) turns into identity 
and the inequalities 1)18(1 are satisfied if 7, z vary in the 



intervals < 7 < 00, — 1 < z < 1. So, the kernel V is 
now parametrized as V = V(7, 2; 7', z'), where 7 is in the 
same interval as 7' and similarly for z, z' . Note that 7, z 
are related to the standard LF variables k±, x used in sect. 
Has j = k 2 ± , z = I - 2x. 
In this way we get: 



7T(1 - Z 2 ) ,/_ 

d/3' 



7 fc 



LU ■ p' 



1-z 2 



+ K 2 -/3'- ie) ( 7 |f z T 2 + k 2 + /?' - ie 



(20) 



The polynomial A(k,p) determining 7(fc,p) has the prop- 
erty: 



A fc 



+ — ,p) =A(k,p)+cp,p', 

LJ-p J 



where cpi = d 2/4(1 — z 2 ) and d is given in appendix 1X1 
Therefore the integrand in 1(20(1 has two poles of the first 
order at (3' — (3± with 



a , ■ 7 + z 2 m 2 2 
P+ = -, 9 K , P- 



7 + z 2 m 2 o 
ie + ^ h K 



1-z 2 ' ' " • 1-z 

and one pole of the third order due to the factor 

1 



lA(k,p)+cp,(3> + ief 

If cri > (i.e., d > 0), the third order pole is at (3' ~ — ie. 
We close the contour in the upper half-plane and take 
the residue at = 0+ ~ +ie. If d < 0, this pole is at 
j3' = +ie. Then we close the contour in the lower half- 
plane and take the residue at j3' = /3_ ~ — ie. The result 
reads: 

W ci )(7,z; 7 ',z') 



7 + m 2 z 2 + k 2 (1 — z 2 ) 

f7(fc + T ^ !?5 ),if C '>0 
\7(fc + ^, J3 ),ifc'<0 



Substituting here expression 1(16(1 for 7, we finally find: 



^ Ci >(7,z; 7 ',z') = -' 



(1-z 2 ) 



2\3 



7T 2 7 + m 2 z 2 + K 2 (l — Z 2 ) 

1 pi rl-V3 /• 1-2/2-1/3 

2 



x / 2/4(1 - 2/4) / dy 3 / 
■/o Jo Jo 



11 



7/ is defined by ((17(1 and 7? reads: 

D = c 7 7 + cy7' + c k k 2 + c m m 2 + c^p, 2 
- y 4 \d\ 7 + m 2 z 2 + k 2 (1 - z 2 ) 



d yi^3' 

(21) 



(22) 



The coefficients determining 7) are given below in ap- 
pendix 
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A.l The coefficients determining D, eq. 1)22(1 

c = 2/2 2 (1 + (-1 + yi + Va) Vi) (1 + z) 

+ 2/3 (Vl 2 2/4 ("I + z) - (-1 + t/ 3 ) (-1 + 2/4) (z - Z') 

+ 2/1 (1 + 2/3 2/4 (-1 + z) - Va z + (-1 + 1/4) «')) 

+ 2/2 ((-1 + 2/4) (1 + z) 

+ 2/i 2 2/4 (1 + 2) + 2/3 2 2/4 (1 + z) 

+ y 3 (1 + 2 2 - z' + y 4 (-2 - 3 z + z')) 

+ 2/1 (l-z'+y 4 (-2-z + 2y 3 (1 + z) + z'))) 

c 7 = 2/4 (2/2 2 (1 + (-1 + 2/1 + 2/3) 2/4) (1 + ^) 
+ 2/3 (-1 + 2/3 + (1 + (-1 + 2/l) 2/l) 2/4 

- 2/3 2/4 + 2/1 2/3 2/4 + 2/1 ^ - 2/i 2 2/4 z - 2/1 2/3 2/4 z 
+ (-l+yi+ 2/3) (-1 + 2/4) 22') 

+ 2/2 (-1 + (1 + (-1 + 2/1) 2/1) 2/4 - z 
+ 2/ 3 2 2/4 (1 + z) 

+ y 3 (2 + z + y A (-3 - 2 z + 2 ?/i (1 + z)) 
+ (-I + 2/4) z^') 

+ z (j/4 + yi (1 - z' + 2/4 (-2 + j/i + «'))))) 

<V = (- 1 + »*) (-1 + (1 + 2/i 2 + (-1 + 2/3) 2/3 
+ 2/1 (-1 + 2J3)) y 4 ) (-1 + z 2 ) 

Ck = - ((-1 + z 2 ) (-1 + z' 2 + 2/4 (2 + (-1 + 2/3) 2/3 

+ 2/2 (-1 + 2/2 + (-1 + 2/2 + 2/3) z) 

- 2 Vl 2 (-l + z')-2/i (-2 + ?y 2 (2 + z) 
+ y 3 (2 + z)) (-1 + z') 

- 2/3 (2^ + (-1 + 2,2 +2,3) ^) z'-2z' 2 ) 

+ 2/ 4 2 (-1 + 2/2 + 2/3 - 2/3 2 + 2/2 2 (-1 + 2/3) (1 + z) 

+ 2/2 (-I + 2/3) (2/3 + (-1 + 2/3) z) 

+ (-1 + 2/3) 2/3"' + 2/2 2/3 (2 + z) z' + z' 2 

+ 2/1 2 (-2 + 2/2 + 2/3 + 2/2 z - 2/3 z + 2 z') 

+ 2/1 (2 + 2/2 2 (l + z)-2z' 

+ 2/3 (-3 + 2/3 - 2/3 z + (2 + z) z') 

+ 2/2 (-3-2z + 22/3 (l + z) + (2 + z) z'))))) 

Cm = -(2/l 3 2/4 2 (-1 + Z 2 )) 

+ 2/2 2 2/4 (1 + (-1 + 2/3) 2/4) (1 + z) (-1 + z + z 2 ) 
+ 2/1 2 2/4 ((1 + z) (-1 + z + 2/2 2/4 z 2 ) 
+ 2/3 (2/4 - 2/4 z 3 ) + 2 (-1 + 2/4) (-1 + z 2 ) z') 
+ 2/2 2/4 (z 2 (-1 + 2/4 + 2/3 (2 + (-3 + 2/3) 2/4) 

- z+ (y 3 + (-1 + 2/ 3 ) 2 2/4) z) 

+ 2/3 (-I + 2/4) (-2 + z 2 (2 + z)) z') 
+ (-I + 2/4) (-((-I + 2/3) 2/3 2/4 z 2 ) 



+ ("I + 2/3) 2/3 2/4 z 3 z' + (-1 + 2/4) (-1 + z 2 ) z' 2 ) 

+ 2/i 2/4 (- (2/3 (2/4 (1 + 2/3 (-1 + z)) - z) z 2 ) 

+ 2/2 2 2/4 (1 + z) (-1 + z + z 2 ) 

+ (-1 + 2/4) (2 - 2 z 2 + 2/3 (-2 + z 2 (2 + z))) 2' 

+ 2,2 (-((-2 + z 2 (2 + z)) (-1 + z')) 

+ 2/4 (2-2 z' + z 2 (-3-2z + 22/3 (1 + z) 

+ (2 + z)*')))) 

c M = (-1 + 2/1 + 2/2) 2/4 (-1 + (1 + 2/i 2 + (-1 + 2/3) 2/3 
+ 2/1 (-I + 22/3)) 2/4) (-1 + z 2 ) 



B Calculating LFD cross- ladders 

At first we will transform the 3-dim. LF equation (jSJ) in 
a 2-dim. form. The kernel Vlf (k'±, x'; k±,x, M 2 ) depends 
on the scalar product 

k± ■ k' j_ = k±k' ± cos (/>'. 

Since t/j(k'±,x) depends on |fe'j_| = k'x and does not de- 
pend on the angle <f>' , one can integrate the kernel over </>' 
and the equation JHJ turns into: 



M ip(k±,x) 



k\ + rn 2 
x(l — x) 

— ( Mh'. ^'\V r „(h'. J h , -r M ^ *'•'"•''< <lr ' 
2tt 3 



(23) 



2.T'(l-a;') , 



where 



2tt 



V LF (A;l,a;',fcj.,a:,M^ = / V LF (k\, x' , k ± , x, M^dfi . 

Jo 

(24) 

To calculate the kernel in the LF equation (JSJ, one can 
use the Weinberg rules [2], which are equivalent to the 
graph techniques in LFD |10|. To calculate the amplitude 
— AA, one should put in correspondence: to every vertex - 
the factor g, to every intermediate state - the factor, 

2 , V- fcu. + m ? 
, where s mt = 



«o - s int + iO 



and so is the initial (=final) state energy. In our case 
(bound state) : Sq = M 2 . To every internal line one should 

put in correspondence the factor e ^ x *^ . One should take 
into account the conservation laws for kn_ and xi in any 
vertex and integrate over all independent variables with 
the measure d f%£y> Xi ■ 

Applied to the ladder graphs fig. |3 these rules result 
in eq. © for the ladder contribution. The integral 124|) is 
calculated analytically. 

There are six cross-ladder time-ordered diagrams shown 
in fig. [3] Sum of them determines, by eq. (|10|) . the kernel 

Vlf ■ Consider four diagrams for Vi_4 in fig. El They 



J. Carbonell, V.A. Karmanov: Cross-ladder effects in Bethe-Salpeter and Light-Front equations 



9 



contain eight lines (a' , a" , a, 6' ,b" ,6, c, d) and three inter- 
mediate states. The momenta corresponding to these lines 
are the following: 

(a) fe X; x 
(a") fc j_, £E 

(a) fex, •£ 
(6') -fc'x, 

(6") fc"x - k' x -fex, 1 + a:" -x'-x 

(b) — fex j 1 — x 

(c) fe'x-fe"x, x'-x" 

(d) fex -fe"x, x-x" 

Introduce for every line the energy: 



E„, 



En 



k' 2 ± + m 2 
x' 

k"\ + ro 2 

x 7 ' 
k 2 ± + m 2 

k'l+m 2 
1-x' 



(k" ± - k\ - fex) 

1 + x" — x' — x 
k\ + m 2 



2 i ™2 



1-x 
(k'x - k"xf + /i 2 

x' — x" 
(fcx-fe"x) 2 +M 2 



Ey — 

E b n = 

E b = 

E c = 

E d = 

x — x" 

Energies in three intermediate states for Vi are written as: 

si = E a „+E c + E b ,, 

s 2 = E a n + E c + E d + E b „, 

s 3 = E a „+E d + E b . 

The kernel Vlf is related to the amplitude M. as: 
Vlf = — -M/(4m 2 ). For the diagram fig. EJVi it has the 
form: 

f g 4 6(x")6{x - x")9(x' - x")0(l + x" - x - x') 
1 ~ J 4m 2 (si - M 2 )(s 2 - M 2 )(s 3 - M 2 ) 

1 d 2 k'ldx" 



x"(x - x")(x' - x")(l + x" - x - x') 2(2tt) 3 



(25) 



d 2 k'ldx" = k'xdk'xd<t>"dx". Therefore the contribution of 
fig- 0^1 to the kernel Vlf, eq. (JUJ, is determined by a 
4-dim. integral: 

V 1 (k'x,x',kx ) x ) M 2 )= / Vi(k'x,x',kx,x,M 2 )d(f)' 
Jo 



i{x,x ) 



dx" 



-4m 2 a 2 k' ± dk'li 



Since: 



max(0,x+x'-l) J n ( s l ~ M 2 )(s 2 — M 2 )(s 3 — M 2 ) 

1 

X x"{x - x")(x' - x")(l + x" - X - X') ' ' 

with a = g 2 / (IQirm 2 ). We have removed the theta-functi- 
ons and incorporated the integration limits in the variable 
x" explicitly. One can calculate the 4-dim. integral (|26|) for 
the kernel numerically. 

Consider now another three kernels V%, V3, V4 shown in 
figs. El They still have the form but the corresponding 
energies s±, S2, S3 are different. 

We denote the kernel corresponding to fig. El V% , after 
integration over <f> ' , as V2(k' ± , x' , fex, x, M 2 ). It has exactly 
the same form as eq. (I2fil) . but with energies given by: 

si = E a i + E c i + E b ", 

s 2 = E a „ +E c + E d +E h „, 

S3 = E a + E c + E b „. 

The kernel V3(k' ± ,x' , kx,x, M 2 ) corresponding to the 
graph fig.0 V3 has exactly the same form as eq. ()26fl. but 
with energies given by: 

si = E a n + E c + E b i , 

s 2 = E a „ +E c + E d + E v >, 

s 3 = E a + E c + E b „. 

The kernel t4(fej_, x ', fex, x, M 2 ) corresponding to the 
graph fig. EI V± also has exactly the same form as eq. (|26|) . 
but with energies given by: 

si = E a i + Ed + E b ", 

s 2 = E a „ +E c + E d + E b// , 

s 3 = E a „+E d + E b . 

Another two cross-ladder graphs are shown in figs. [31 
V5 and Vq. Here we have new lines d',c'. Corresponding 
momenta and energies are the following: 

(d!) k" x fex, x" -x 
(c 1 ) k"x - k'x, x" - x' 



(k"x-kx) 2 +^ 2 

&d' — 7, — 

x" — x 
(k"x-k'x) 2 + » 2 

I^r' — 7. : = — £Jr 



fex • fe'x = kxk' ± COS(f>', 

fex ■ fc"x = kxk'l cos</>", 
k'x k"x = k'xk'i cos(^'-0"), 

the integrand depends on two azimuthal angles <p',<f>". 
The kernel 1)25(1 is determined by the 3-dim. integral over 



If x' > x, only V5 contributes. It has the form: 

-4m 2 a 2 d 2 k'ldfi 



Vr = / dx 



7t(«i - M 2 ){s 2 - M 2 ){s 3 - M 2 ) 
9(x' - x) 



x"(x" - x)(x' - x")(l +x" - x - x') 



(27) 
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with si, S2, S3 given by: 

si = E a »+E c + E v , 

52 = E a + E c + Ed' + Ey , 

53 = E a + E c + Eyi . 

If x > x', only Vq contributes. It has the form: 

-4m 2 a 2 d 2 k"dd>' 



V 6 



dx ' 



tt{ Si -M 2 )(s 2 -M 2 ){s 3 -M 2 ) 
6{x-x') 



(28) 



x"(x - x"){x" - x')(l + x" - x - x 1 ) 

with si, S2, S3 given by: 

si = E a < + Ed + Eb" , 

s 2 = E a ^ + E d + E d + E b , 

s 3 = E a „ +E d + E b . 

B.l Test by the Feynman cross graph 

On the energy shell the sum of six contributions I/( CL ) — 
Si=i ^ %■ El must coincide with the value of the Feyn- 
man cross graph, fig. ^ taken on the mass shell, namely 
with 

i^ (Ci) (M>). 

K^ CL \k, k' ,p) is given by (|13[) where in the denominator 
A' , eq. I|14|) . one should put 



M 2 = s, k 2 = k' = m 2 - -s, 

k ■ p = k' ■ p = 0, k ■ k' = m 2 s 1 

4 2 

and s, t are in the physical domain. Expression i|14|) is 
considerably simplified and we find: 



V, 



(CL) 



Feyn 



-4a 2 m 2 



1 /-i~y3 pi— V2— vs 

dy 3 / dy 2 / dyt 
Jo Jo 



(i-yi- 2/2) - y3t) + (2/1 - 2/2) to 



+2/3* + 2/12/2 (s + i) 



(29) 



The 3-dim. integral l|29l) is calculated numerically. One 
can chose a particular kinematics, for example, in c.m.- 
frame, where s and t for elastic scattering are determined 
by incident (=final) momentum and the scattering angle. 

On the other hand, we calculate the sum of six LF 
cross graphs 



V, 



LF 



{k' x ,x',k x ,x,M 2 = s ) = £Vi 



(not integrated over 4>'). Here 

^ _ uj - pi 

ui ■ (pi + p 2 ) ' 



uj-p 1 



w ■ (Pi +P2) 



and k±,k'± are projections of pi, p'i on the plane or- 
thogonal to ui. We remind that uj — (ct>o, <jj) is an arbitrary 
four-vector with uj 2 = 0. Though, for a given kinematics, 
the values of k±, x, k'±, x' depend on orientation of uj, 
the on-shell amplitude vj^, L \k' x' , k±,x, M 2 ~ s) does 
not depend on it. We can chose us = (1, 0, 0, —1), find the 

arguments of vj^ L \k' j_, x' , kx, x, M 2 — s) in kinematics 

for the cm. elastic scattering and calculate vj^ L \ 

For different incident momenta and scattering angles 

we have checked numerically that Vp^^, eq. I|29|) . coin- 
cides with the sum of six cross graphs vj^ L ^ — J2i=i ^5- 
This confirms validity of both calculations. 



C Calculating LFD stretched boxes 

The stretched-box graphs determine the kernel V£ F , eq. 
(|10|l . They are shown in fig. p|| and their contributions are 
denoted as V7 and Vs- We introduce the line b'" '. All the 
momenta except for the line b'" are the same as for the 
cross graphs. The momentum and energy carried by the 
lines b'" are the following: 

(6"') -k" x ,l-x" 



„//2 



Eyn 



m 



1-x" 



The following product of the theta- functions enters the 
kernel V7 

6{x")9{x" - x)0(x' - x")6{l - x"). 

It restricts the integration domain as: 

< x < x" < x' < 1. 

Therefore the kernel V? is nonzero at x 1 > x only. It has 
the form: 



V 7 



dx" 



-4m 2 a 2 d 2 k''d<j)' 



7r(s! - Af 2 )(s 2 - M 2 )(s 3 - M 2 ) 
9(x' - x) 



(30) 



x"{x" - x){x' - x")(l - x") 

The corresponding energies s read: 

si = E a „ +E c + E b >, 

s 2 = E a + E d , +E C + E V , 

s 3 = E a + E d , + E h n, 

and the energies E are given above in appendix and in 
this appendix. 

The kernel Vs is non-zero at x' < x only. It has the 
form: 



^8 = 



dx" 



-Am 2 a 2 d 2 k"dd>' 



n(si - Af 2 )(s 2 - M 2 )(s 3 - M 2 ) 
6{x-x') 



x"[x-x"){x" -x')(l-x") 



(31) 



J. Carbonell, V.A. Karmanov: Cross-ladder effects in Bethe-Saipeter and Light-Front equations 



with si, S2, S3 given by: 

si = E a ,+E c ,+E v „, 
s 2 = E a > +E d +E d + E b , 
s 3 = E a n +E d + E b . 
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